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Abstract 


Mass transfer of gaseous components from rising bubbles to the ambient 
liquid can be described based on continuum mechanical sharp-interface 
balances of mass, momentum and species mass. In this context, the stan¬ 
dard model consists of the two-phase Navier-Stokes equations for incom¬ 
pressible fluids with constant surface tension, complemented by reaction- 
advection-diffusion equations for all constituents, employing Pick’s law. 

This standard model is inconsistent with the continuity equation, the mo¬ 
mentum balance and the second law of thermodynamics. The present 
paper reports on the details of these severe shortcomings and provides 
thermodynamically consistent model extensions which are required to cap¬ 
ture various phenomena which occur due to the multi-physics of interfacial 
mass transfer. In particular, we provide a simple derivation of the interface 
Maxwell-Stefan equations which does not require a time scale separation, 
while the main contribution is to show how interface concentrations and 
interface chemical potentials mediate the influence on mass transfer of a 
transfer component exerted by the change in interface energy due to an 
adsorbing surfactant. 

Keywords: Soluble surfactant, interface chemical potentials, interfacial en¬ 
tropy production, interface Maxwell-Stefan equations, jump conditions, surface 
tension effects, high Schmidt number problem, artificial boundary conditions. 

Introduction 

The standard model for the continuum mechanical description of mass transfer 
across fluid interfaces is based on the incompressible two-phase Navier-Stokes 
equations for fluid systems without phase change. To be more precise, inside 
the fluid phases the governing equations are 


V-v = 0, 

9t(pv) -I- V ■ (pv (g) v) -I- Vp = V • -t- pg 


IVISC 


( 1 ) 

( 2 ) 


with the viscous stress tensor 


S 


I Vise 
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where the material parameters depend on the respective phase. Whenever dis¬ 
tinction between the different phases is necessary, -|- and — are used as phase 
indices. Since phase change is neglected in the standard model, there are no con¬ 
vective fluxes across the interface, i.e. the normal component Vs = -riE of the 
interfacial velocity coincides with the normal component of both the adjacent 
fluid velocities. Moreover, since both fluids have non-vanishing viscosity, no-slip 
between the two phases at the interface is usually imposed. Finally, constant 
surface tension is assumed. In this situation, the interfacial jump conditions for 
total mass and momentum are 

M = 0, |-S] ■ ns = o-Ksns, (4) 

where S = —pi -f is the stress tensor and ks = —V • ns is twice the mean 
curvature of the interface. Here ns is the unit normal at the interface directed 
into the bulk phase -, say, and the notation 

|0](x) = lim (0(x-f/ins) - ^f)(x -/ins)) (5) 

stands for the jump of a field (j) across the interface. 

The local molar concentration Ci of a chemical species Ai is governed by the 
balance equation 

dtCi + V ■{ci-v + 3i) = ri, ( 6 ) 

where the molecular fluxes are typically modeled according to Pick’s law as 

Ji = -DiVci (7) 

with constant diffusivity. The source term on the right-hand side in ([S]) accounts 
for chemical reactions. At the interface, the diffusive fluxes in normal direction 
are commonly supposed to be continuous, i.e. 

l-AVc,] • ns = 0. (8) 

One more constitutive equation is needed to determine the concentration profiles 
of Ai, where instantaneous local chemical equilibrium at the interface is usually 
employed. This means continuity of the chemical potentials at the interface, i.e. 

b.] = 0. (9) 

Using standard relations for the pi, this leads to Henry’s law. In the form for 
molar concentrations, the latter states that 

Cfe = ct/H, (10) 

with a Henry coefficient Hk, often assumed to be constant. Equations o- 
and (l51)- (fTU)) comprise what is called the “standard model” throughout this 
paper. This standard model, sometimes with further simplifications like homo¬ 
geneous gas phase concentrations or even constant liquid-sided concentration at 
the interface, is the basis of almost all detailed numerical simulations of mass 
transfer across fluid interfaces up to now; see m, [m and the extensive list of 
references given there. For fundamentals on continuum mechanical modeling of 
two-phase fluid systems we refer in particular to |2S1 [T71 [21 [32] . 
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The remainder of the present paper reports on major short-comings of the 
standard model and explains how to overcome these in the framework of contin¬ 
uum thermodynamics. For this purpose, we start on the interface, proceed to 
the adjacent bulk layers before we visit the bulk phases and, finally, the outer 
boundary of the domain in which the balance equations are studied. To enable a 
more descriptive notion, we focus on mass transfer from a bubble to the ambient 
liquid as a prototype case, but all considerations are also valid for other phase 
topologies (drops, films etc.) and for liquid-liquid systems. 


The Fluid Interface 

In our continuum thermodynamical description of the interface, a sharp-interface 
model is employed which represents a thin layer in which the partial mass den¬ 
sities change from one to the other (local) bulk value within a layer of thickness 
S in the order of a few Angstrom. Since every transfer component of the 
mixture, in principle every constituent, is also present in this layer, it has a 
non-zero interface concentration cf, hence also an interface chemical potential 
Ilf in the sharp-interface model. The mass balance for Ai on the interface then 
becomes more involved and reads 


cf -f Ve • (cf -k Jf) -k [cAv - v^) -k J*] • riE = rf, (11) 

where denotes the interfaces barycentric velocity, Jf the interfacial diffusive 
flux and rf is the total molar rate of change of Ai due to interface chemical 
reactions between the species. Furthermore, dfcf denotes the time derivative 
- here of interface concentrations - along a path which follows the interface’s 
normal motion; note that a partial time derivative of interface quantities does 
not exist in the usual sense if the interface moves. Observe that the surface 
divergence in m contains a curvature dependent contribution since 

Ve ■ = Ve ■ - bsKEns, (12) 

where the subscript || denotes the tangential part. Note also that even for 
vanishing interface concentrations and without interface chemistry, equation 
m does not reduce to (|5]) from the standard model, but the resulting jump 
condition also contains a convective part due to a relative motion of the interface 
to the bulk matter. This is crucial, for example, to model the condensation of 
a vapor bubble or the dissolution of a pure gas bubble. 

If (ITTl) is multiplied by the molar mass then the partial mass balance 

9fpf + Ve • (pfv^ -k jf) -f IpAv - v^) -k ji] • He = M^rf (13) 

results with jf = M^Jf the diffusive mass flux. It is helpful to write out the 
jump-bracket, i.e. 

IpAv-v^)-kji]-nE = -(p+(v+-v^)-kj+)-n+-(p-(v--v^)-kjA)-n- (14) 


with denoting the outer unit normals to the respective bulk phases. We 
abbreviate the terms on the right-hand side by letting 


m. 


±,s 




(15) 
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since these terms denote the mass transfer rates from the respective bulk phase 
to the interface and comprise the central object of the present paper. Guided 
by the case of a (soluble) surfactant, i.e. a chemical species which significantly 
accumulates at the interface, thereby changing the interface energy, the one¬ 
sided mass transfer terms are further split into an adsorption and a desorption 
term according to 




ad,-\- de,-\- 


• —,E ad,— 
m- =S, 


de,— 

— s, . 


The interfacial balance for total mass reads as 


(16) 


dfp^ + Ve • (p^v^) -f |p(v - v^)] • riE = 0, 


(17) 


and we abbreviate the total mass transfer rate from the respective bulk phase 
to the interface as = p^(v^ — v^) • n^, while inside the jump bracket m 

stands for m = p(v — v^) • ns with corresponding one-sided limits, but with ns 
instead of n^. Summation of all partial mass balances for the interface hence 
shows that 

N 

= o ( 18 ) 

i=l 

always holds for the diffusive mass fluxes. The latter implies a corresponding 
constraint for the diffusive molar mass fluxes which is not satisfied by Fickean 
fluxes according to ©• 

Equation (I17|l replaces the normal component of the first equation in ([5]), 
while the general form of the second equation in (jT]) is the interfacial momentum 
balance 


df ) + Ve • (p^ ® ) + |pv® (v - ) - SI • ns = Ve ■ (19) 

with specific body force and the surface stress tensor S^ which contains the 
surface tension cr, commonly used instead of the surface pressure = — cr, and 
viscous interface stresses due to intrinsic surface viscosities. For a fluid interface 
with Newtonian rheology this leads to the Boussinesq-Scriven surface stress, i.e. 

S'" = (a + (pf - pf )Ve • v^)Pe + 2pf(20) 

with the interface dilatational and shear viscosities pf > pf > 0, the surface 
projector 

Pe = I — He G He (21) 

and the surface rate-of-deformation tensor = ^Pe(Vev^ -|- (Vev^)^)Pe- 
For further information about interfacial rheology we refer to while m the 
present work we rest content with a reduced form which assumes that, while the 
interface contributions in the partial mass balances are of crucial importance, 
their inertia can be neglected and that no intrinsic interface viscosities appear. 
Then the interfacial momentum balance becomes a momentum transmission 
condition which reads as 


|rhv] = |S] • riE + o-KEns -b Vect. 


( 22 ) 
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Note that (I22|) is an interfacial force balance which accounts for the effects from 
mass transfer (the Stefan flow on the left-hand side) and for Marangoni effects 
expressed by the last term on the right-hand side. If the fluids are stagnant, i.e. 
if no flow occurs, then (EH) reduces to 

bJns = ctkeiie -f Vect. (23) 

The normal component of (l23ll implies the Young-Laplace law, i.e. 

p+ -p- = ‘^ (24) 

in case of spherical bubbles or droplets of radius R, while the tangential part im¬ 
plies constant surface tension (individually on connected interface components), 
since there is no term to balance the tangentially acting Marangoni stress. The 
important consequence of this is 

Vect ^0 V 7^ 0, (25) 


i.e. any Marangoni stress induces a flow at the interface, the so-called Marangoni 
convection, which usually has a significant impact on mass transfer processes. 

To obtain a closed model, i.e. a system of partial differential equations con¬ 
taining rates and fluxes which are expressed in terms of the primitive (balanced) 
quantities up to some material parameters which must be determined from mea¬ 
surements or a micro-theory, one needs constitutive equations for the unknown 
rates and fluxes. At the interface, these include the individual one-sided mass 
fluxes rhj ’ . Moreover, one needs constitutive equations for the interface chem¬ 
ical potentials and an interface equation of state modeling the surface tension. 
The latter quantities all follow from a constitutive equation for the interface free 
energy via the surface Gibbs-Duhem equation 

N 

=a + Y^ pfpf (26) 

and the fundamental differential relation 

d{p^'ip^) = -p^s^dT^ -I- ^ dfdpf, (27) 

i=l 


where denotes the specihc interfacial entropy. The latter implies in particular 

E _ d{p^^^) 


= 


dpi 


for p -tp = p i!) {T ,pi,...,p^). 


(28) 


The constitutive equations should of course be consistent with the second law 
of thermodynamics, i.e. such that the entropy inequality is satisfied: ^ 0 

for any thermodynamic process. Now, in the considered case with partial mass 
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densities on the interface, the interfacial entropy production is given by 



A derivation of (1^^ along the lines of classical TIP can be found in [5] , where the 
relation (EH) is imposed. The interfacial entropy production has been derived 
in different form before in m, based on the single component case treated 
in [3]. The Gibbs-Duhem relation (l26ll can also be obtained as a consequence 
and we briefly explain the main steps in an appendix. The interfacial balances 
for partial mass, total momentum and energy can also be found in |32| and 
El- While in m the entropy production is only considered in the special case 
with pf = 0, the interfacial entropy production in m is differently structured 
and the one-sided partial mass transfer rates rh^ ’ do not appear as separate 
quantities. A similar form of the interfacial entropy production can be found 
in EI1I3II, again without terms containing the one-sided partial mass transfer 
rates At this point it is also important to note that in im [m [311132] 

the concept of Gibbs excess quantities is employed, which is known to lead to 
problems like negative excess of interface partial mass. Employing this concept, 
a quantity which is not accumulated at the interface would typically receive no 
excess mass which implies that its mass transfer will not involve an interface 
chemical potential. But the latter is shown below to account for the influence 
of changes in the interfacial free energy. 

Each binary product in (1291) corresponds to a dissipative physical mechanism 
at the interface. In each product, one factor is a constitutive quantity which 
is to be modeled in dependence of its co-factor. For processes not too far from 
equilibrium, a linear (in the co-factors) closure is appropriate. In this way the 
first product, which relates to heat conduction, leads to Fourier’s law for the 
interfacial heat flux, i.e. 

qS = -A^VeT^ with A^ > 0. (30) 

The interface temperature is, in general, different from the one-sided limits 
of the bulk temperature. Note also that A^ is allowed to depend on and 
the pf. The latter is also true for all further closure parameters which are 
introduced below without explicit mentioning. 

The second product in (1291) corresponds to multicomponent diffusion on the 
interface. The standard linear closure requires first to incorporate the constraint 
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m, say by eliminating the flux j^. This leads to the mass diffusion fluxes 



where the matrix [Lik\ of mobilities is positive definite. It can also be shown that 
{Lik\ is symmetric, i.e. the Onsager symmetries hold; see [HI for a rigorous theory 
of Onsager symmetries for transport coefficients. A difficulty with this general¬ 
ized Fickean form is the fact that the mobilities are complicated functions of the 
composition - this needs to be the case, since for constant mobilities the posi¬ 
tivity of solutions would not be guaranteed. For this reason the Maxwell-Stefan 
form of surface diffusion is advantageous, and we provide a simple derivation 
based on (I?7ll and being fully rigorous in the isothermal case. At this 
point it should be noted that the diffusion flux is related to a diffusion velocity 
uf via 

jf = pfuf ■= pfi^f - (32) 

where is the continuum mechanical velocity of the individual component Ai. 
Then the starting point is the following reformulation of the relevant entropy 
production contribution. We have 


N 


c, 


DIFF 


= - II 


2=1 


„EV7 _ Pi A 

Pi rp^ rpY: Pi ^ 


(33) 


with a Lagrange parameter A which can be arbitrary due to (I18II . This param¬ 
eter is chosen in such a way that the co-factors in (|33|l sum up to zero, hence 



Using (1^51) and (ITfl) . a straightforward computation shows that (155)) then be¬ 
comes 


N 


Cdiff = - H uf • 


(35) 


with the generalized thermodynamic interface driving forces 

'11^ 


di Pi Pi 




(36) 


where yf = pf /(P are the interfacial mass fractions, 2 /fbf the total 

interface body force density and YiP is the interface enthalpy density. From a 
kinetic theory, or from a more elaborate continuum theory employing partial 
momenta (see 0), the interface driving forces are 


jE _ th _ U 

di Pi Pi rpT. 


^ h^\7 ^ 


(37) 


with the partial interface enthalpy densities hf. Note that both agree in the 
isothermal case. Now, equation (I35p is used to derive constitutive equations for 
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the df, after , say, is eliminated by means of the constraint 


N 


Edf = 0. 


( 38 ) 


Linear closure yields 

A^-l 

df = - ^ rifc(u^ - u^) for i = -1 (39) 

k^l 

with a positive definite matrix [Tik] of interaction coefficients. By arguments 
fully analogous to those given in [B] for the bulk diffusion case, the assumption 
of binary interactions, i.e. Tik = Tik{T^, pf, Pk) 0 if P? Pfc to 0, leads 
to the Maxwell-Stefan form 

N 

df =-'^fikpfpkiuf-uf) for i = (40) 

fe=i 

with positive, symmetric friction coefficients fik- An equivalent form for molar 
quantities is 


-E 


tS 


Jf- 




Dfk 


= df 


for i = 1,..., N 


(41) 


with the interfacial Maxwell-Stefan diffusivities It turns out that inversion 
of the system m together with (fTHll leads to diffusion fluxes which guarantee 
positivity of the solutions of the final partial differential equations for the species 
concentrations; cf. [3 [21]. The interface Maxwell-Stefan equations are also 
derived in but under the assumption of quasi-stationary hydrodynamics 
which is not needed here. 

The third term in (PI represents the entropy production due to interface 
chemical reactions, where iij is the rate of the a*^ reaction and the associ¬ 
ated affinity. For chemical reactions, a linear closure is usually not appropriate. 
Since interface chemistry is not in the focus of the present work, we refer to 
|B| for a non-linear closure of chemical reaction rates which, in simplest cases, 
resembles the well-known mass-action kinetics. 

The 4:^^ and S**' binary products in correspond to the dissipation due to 
momentum transfer between the bulk phases and the interface. Linear closure 
yields the Navier slip conditions 


(v^ — v^)||-I- q;"'^(S^ . n^)|| = 0 with a"*" > 0. (42) 

A common choice is = 0 which leads to continuous tangential velocities at 
S as used in in the standard model. 

The 6*^ and 7*** binary products in (l2^ refer to entropy production during 
energy transmission from the bulk phases to the interface. Linear closure yields 
the relations 


1 1 
j'± 



with /3^ > 0. 


(43) 



A common choice is /3^ = 0 which leads to continuous temperatures at S. 

The final two binary products in (1291) describe the entropy production due to 
transfer of partial mass across the interface. If cross-effects between the species 
are ignored, linear closure gives 


±.E 


= 7i 


hL 

T± 


yS 


1 

J'T, 


(v±-v^) 


i :\2 


— n 


g±,visc 


n 


(44) 


with 7 )^ > 0. Equation shows that mass transfer is driven by chemical po¬ 
tential differences but, in addition, a kinetic term and viscous forces contribute 
to the driving force. In the limiting cases, as > 00 , one obtains 


Pf Pf 1 ^ 

T± y 2 


(45) 


These limiting cases are referred to as those of vanishing (one-sided) iiiterfacial 
resistance against mass transfer. If the temperature is assumed to be continuous, 
often an appropriate assumption, and if the kinetic and viscous terms can be 
neglected compared to the chemical potentials, then (1451) yields 

pt=Pi=pf'^ (46) 


cf. [5] for an estimation of the strength of the different contributions in (H51) . 
Let us note that the first equation in (HSl) . which is usually used to describe 
the concentration jump at the interface, has to be evaluated with care, since 
not only the respective one-sided limits of bulk compositions enter the chemical 
potentials, but also the different pressures: 

PiiT,p^,xy,.. .,x^) = fi~{T,p~,x ^,... ,x((f). (47) 


Since the pressure also has a jump at S and the height of this jump depends on 
the curvature, the pressure dependence of the chemical potentials introduces a 
curvature influence into (1471) . implying in particular that smaller bubbles display 
an increased solubility. In the context of thermally driven phase transfer, the 
same effect is described by the Kelvin equation. 

Note that if (HH)) or (HTl) is employed instead of (1441) . say, than the mass 
transfer rates are no longer explicitly determined from the driving differences 
in chemical potentials, but only follow indirectly from the bulk concentration 
profiles at the interface for the diffusive contribution and from the total mass 
transfer rate for the convective part. If rh is also not known explicitly, which 
is usually the case since continuity of the temperatures is assumed, then the 
convective part can still be accounted for, but follows from the solution of a 
linear system of equations; cf. muni. 

In order to understand the influence of surface coverage by surfactant to 
the mass transfer of another chemical species, we have to stick to a general 
closure like (HI) . But we prefer a non-linear closure in analogy to chemical 
reactions, viewing the mass transfer from the bulk phases to the interface as ad- 
and desorption processes, since the system can be far away from mass transfer 
equilibrium. For technical simplicity we neglect the kinetic and viscous terms, 
i.e. we exploit the reduced products 

1 ^ 

C4ans = y E ’^) (48) 

2=1 
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to derive the closure relations 


^k:± = ^ - f^f) ^ 0. (49) 

Note that one of the rates, either the ad- or the desorption rate, has to be 
modeled based on a micro-theory or experimental knowledge. Then the other 
rate follows from (03). Below, we will let af = 1 for simplicity which already 
suffices to obtain interesting new results. We first consider the case of a soluble 
surfactant which, for simplicity, is only present in say, and on S. Desorption 
is often more easy to model, where the simplest rate function is 

sf = kf^xf (50) 

with the interfacial molar fraction xf = cf /c^. According to with af = 1, 
the associated adsorption rate is 


ad 


= k° 


de 


xf 


exp( 






RT 


)■ 


(51) 


To achieve a concrete result, assume ideal mixtures both in the bulk and on the 
interface, i.e. 


= gf{T,p) + RT\a.xf (52) 

with g^(T,p) denoting the bulk Gibbs free energy of component Ai under the 
temperature and pressure of the mixture and 


S /ry-i 1 

{T,P 


x^ 

! 1 


, a;^_i) — gf{T,p^) + RTlnxf 


(53) 


with gf{T,p^) denoting the surface Gibbs free energy of component Ai under 
the temperature and surface tension of the mixture. Insertion of the chemical 
potentials into (IHT|) yields 


ad 


k° 


de 


exp 


t 9t - 9 } 

'' RT 


)xt =■ 


pad + 
A/,' O/,- , 


(54) 


where kf‘^ depends in particular on the surface tension. Together, this yields 
the simplest ad- and desorption rates, leading to the so-called Henry isotherm. 

Next, we consider a transfer component Ai, like a dissolving gas, which 
does not accumulate at the interface as a surfactant does, but nevertheless 
needs to pass through the transmission layer at the phase boundary which is 
mathematically represented by the sharp interface. In this case, the interface 
concentration Ai is assumed so small that (IT51) can be approximated by 


Irhij =0 + rh^= 0 . 

Consequently, employing (HU and (|i3 "^ith = 1, we obtain 

,+ 


de,-\- 

S, 


exp ( 


N 


9-f 


RT 


)-l + 


de,— I / 9i 
s, ’ I exp (— 




RT 


)-l =0 


(55) 


(56) 


Employing the same simplifying assumption of ideal mixtures and (|50|l . this 
implies 





9t - 9f 

RT 



+ k 


de,— 



9i -gf 

RT 



0, (57) 
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hence 


xf = 


fc,"'+exp(^^^)x++ exp(^^^)x: 

rfe.+ ^ . de.- 


(58) 


Inserting this value into the first expression in (l57ll . equation dm yields 




r^(ie,+ ide,— 


3 exp (--!-)( exp (f-)x+-exp (f- 


,de,+ 


+ k 


de,- 


RT' 


9i 


9i 


RT^ 


'RT^ 


(59) 


Let us compare (IKUll with a closure which does not account for the interface 
concentrations. In this case, (I55|l is directly build into the entropy production 
such that only a single binary product per species remains, namely 


Ctrans 


1 • IT (v - v^)^ 

~T ^ - 2 - 

i=l 


gvisc 

- ns- 

P 


ns]- 


(60) 


Neglecting again the kinetic and viscous terms, the non-linear closure analogous 
to (HlTll yields 

TO+’^(= -to“’^) = ki exp ( - l^exp {^)xt - exp {^)x~ ^ , (61) 

where the ideal mixture assumption is used. The most important difference to 
(151?)) is that there, the mass transfer rate is influenced by the surface tension via 
the surface Gibbs free energy. In contrast to dSU, this allows to account for the 
effect of surfactants on the mass transfer of the considered transfer component. 
For example, assume the surface equation of state to be given as 

p^ = RTj2cf + K^(4--^)- (62) 


The basis for this equation of state is that the interface with its surface tension 
in the clean state is built by component A]\[ (the solvent, say) as the phase boun¬ 
dary between a liquid and a vapor phase. The phase boundary is modeled as a 
compressible interface phase with compressibility K^. Additional components 
Ai,... ,Am-i are (potentially) present in both bulk phases as well as on the 
interface. Then the surface analog of the construction of a consistent free energy 
as explained in §15 in yields the corresponding interface free energy as 


N-l 


^ref ^ref 


2=1 


2=1 


N 

2=1 


cr In xf^. 


(63) 


The corresponding (molar based) chemical potentials then are 

=RT\n{l + j^)+RTlnxf for i = 1,..., iV - 1 (64) 

and 

Ma = ^ln(l +T^rina;^. (65) 

^ref 
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Hence the ideal mixture form according to (1551) results with 


gf{T,p^) = i?rin (1 + ^) for * = 1,..., iV - 1. 
Insertion of (P|l into ((551) implies the relation 


( 66 ) 


• +.S 
m- = 


ide^+ide,— 1 / "I" ~ 


■ (67) 


Taking the clean surface as the reference state, this yields 

•contain _ clean_ ^ <^clean clean 

i'Ia - II In 


l+P?ontam/7^^ 


- 


^contain 


( 68 ) 


for the mass transfer in a system contaminated by surfactant, given as a relation 
to that in the clean system. The specific relation ((551) evidently results from 
very strong assumptions which are not realistic in particular for high surfactant 
concentrations. If the interface (or parts of it) is covered by a densely packed 
surfactant layer it becomes (locally) almost incompressible. In this case, if the 
interfacial mass density is assumed to be constant, the surface pressure depen¬ 
dence of the surface chemical potentials becomes linear. Then, if the interface 
mixture is still assumed to be ideal, the (molar based) chemical potentials are 
of the form 

M? (T, p^,xf) = fhT) + M,?^+RT\nxf. (69) 

P 

Hence (1551) again holds but with 


gf(r,p^) = /f(r) + M,4- (70) 

p 

Via ((59() . this corresponds to a reduction of the mass transfer rate by an exponen¬ 
tial damping factor of Boltzmann type, i.e. a factor of the form k exp{—ap^/RT), 
in accordance with the energy barrier model due to Langmuir; see m, [la and 
the references given there. Experimental data for the transfer of CO 2 from 
Taylor bubbles under the influence of different surfactants in m also support 
the fact that the surface pressure of the contaminated system - not, in the first 
place, the surfactant concentration - determines the mass transfer reduction. 

The description above provides a means to include mass transfer hindrance 
induced by surface coverage into detailed numerical simulations: if the interfacial 
resistance of the clean system against mass transfer is negligible, say, one can 
compute the local mass transfer rates from the bulk concentration gradients 
at the interface, accounting for Henry’s law to model the concentration jumps, 
and then correct the transfer rates via (1551) or a more elaborate modeling using 
(155() and appropriate interface free energies. Notice that a second significant 
effect of surface coverage is the Marangoni-convection resulting from the 11011 - 
homogeneous surface tension in (1551) which often leads to an immobilization of 
the interface at the rear end of moving bubbles. For spherical bubbles, the latter 
is typically accounted for by the stagnant cap model; for this as well as for a 
more detailed but still simplified model of the immobilization effect see m and 
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the references given there. In future work we will further elaborate on the above 
approach to describe mass transfer under the influence of adsorbed surfactant, 
considering more involved models of the interfacial free energy, also taking into 
account electrical effects due to, say, ionic surfactant or salts. 


The Vicinity of the Interface 


The immediate consequence of the processes at the interface is the occurrence 
of possibly strong gradients in the adjacent bulk regions. This concerns both 
phases and can lead to local depletion of a certain species due to fast transfer 
through the interface into the other bulk phase. In this context it should be 
mentioned that mass transfer problems for multicomponent systems are always 
conjugate problems. This is nicely demonstrated by the experiments on CO 2 - 
transfer from Taylor bubbles into water in m in which case a pure C02-bubble 
is present initially, but during its dissolution into the water, other gas compo¬ 
nents dissolved in the water phase lead to mass transfer into the bubble such 
that it does not completely disappear but a gas bubble, free of CO 2 , finally 
persists. 

The dissolution of a chemical component in the bulk phase into which it is 
transferred typically leads to a release of heat (heat of mixing) which can induce 
locally non-isothermal conditions. This effect is of course stronger for smaller 
bubbles since the specific area scales as I/i? for a spherical bubble of radius 
R. This may explain that small C02-bubbles with diameter below 100 mm 
show faster rise velocities than predicted by the classical Hadamard-Rybczinski 
solution although the latter assumes a fully mobile interface; see m. where a 
local reduction of the liquid viscosity due to a locally increased temperature is 
mentioned as a possible explanation. This effect will be even much stronger in 
case of reactive mass transfer with fast exothermal chemical reaction, since in 
this case the chemical reaction will occur mainly in the region near the interface. 

The presence of fast chemical reactions will also amplify another well-known 
phenomenon which is a severe challenge for so-called direct numerical simula¬ 
tions, i.e. the approximate numerical solution of the mathematical model with 
such a fine resolution in time and space that all relevant time and length scales 
are sufficiently resolved for an accurate and grid independent solution. This 
phenomenon is the occurrence of extremely thin concentration boundary layers 
due to advection-dominated transport. To explain this in some more detail, 
consider the dimensionless form of the bulk transport equation for molar mass. 
Here it suffices to consider the simplest form from the standard model, i.e. 


d;c* + V* ■ ic*v*) = 


I 


ReSci 


(71) 


where chemical reactions are not included for simplicity. The dimensionless fac¬ 
tor in front of the diffusion term is much smaller than 1 in most technically 
relevant cases, showing that advective transport is much stronger than diffu¬ 
sive one. But even more severe, the factor in front of the Laplacian is much 
smaller then 1 /Re, which governs the diffusion of momentum. Recall that typi¬ 
cal Schmidt numbers for dissolved components in liquids are of the order of 1000 
or even a few magnitudes larger for macro-molecular components, e.g. proteins 


13 



or large surfactant molecules. The effect of the Schmidt number is the occur¬ 
rence of finest structures in the concentration fields, like local lamellar sheets 
or boundary layers, with typical length scale ^diff down to the Batchelor length 
scale given as 


^diff 


^conv /v^, 


(72) 


where 4onv denotes the smallest length scale of the velocity field. This was 
derived in [Tj for scalar transport in turbulent flows, in which case Iconv is the 
Kolmogorov length, but is also valid as a lower bound for the length of the 
smallest persisting structures in a scalar distribution as shown computationally 
in m for mixing in a T-shaped micromixer. At present, direct numerical simu¬ 
lations of mass transfer at gas bubbles, say, in fully three-dimensional cases even 
for moderate Reynolds numbers are not feasible for realistic Schmidt numbers. 
This is true even more for the case of reactive mass transfer. The solution to 
this “High Schmidt Number Problem (HSNP)” is one of the largest numerical 
challenges in this field. 


The Gas Phase 

In virtually every two-phase fluid system, the gas phase will be composed of 
more than one chemical component; even if a pure gas bubble is placed in 
a degassed liquid, the gas phase will receive liquid vapor during the bubble 
evolution. Often three or more constituents will be present in the gas phase 
in which case diffusion cross-effects can become relevant. This is a well-known 
fact which leads to phenomena like up-hill diffusion as shown in the classical 
experiment by Duncan and Toor m- A proper description of multicomponent 
diffusion with cross-effects and non-ideality effects is possible via the Maxwell- 
Stefan equations together with a proper modeling of the chemical potentials 
including especially their pressure dependence. The bulk version of the Maxwell- 
Stefan equations are completely analogous to dT51) . (1571) and (ITT]) above, only 
that (7 is to be replaced by —p in (EZD. For more details on the Maxwell-Stefan 
description of multicomponent diffusion see [SI31I1]. A second important weak 
point of the standard model is the assumption of isochoric flow inside the gas 
phase, i.e. zero divergence of the barycentric velocity field. This excludes volume 
variations of the bubble as a reaction to changes in pressure, but the latter can 
be significant in particular for bubbles rising due to buoyancy, since such bubbles 
can experience strong variations of hydrostatic pressure which typically leads 
to expansion of the bubbles during their rise, accompanied by a change of the 
constituents’ chemical concentrations. To include such pressure effects into the 
modeling, one should note that incompressibility of the gas phase is usually 
introduced via the low Mach number (Ma) approximation. The latter employs 
an expansion of the balanced quantities as power series in Ma and, by equating 
coefficients for corresponding powers, the pressure contributions pq and pi are 
ruled out, while p 2 is the pressure which remains as a Lagrange parameter in 
m, taking care of the zero divergence condition. But a closer look at the low 
Mach number approximation shows that po and pi are not constant: They have 
zero gradients, but are still allowed to depend on time. This indicates a possible 
way to include the pressure dependence inside the gas phase. 
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The Liquid Phase 

As in the gas phase, the assumption of incompressibility in the standard model 
can also be a weak point for the liquid phase. First, incompressibility is often 
even imposed via the assumption of constant mass density. At this point, note 
that the physical definition of incompressibility for a single component material 
states that the mass density given as a function of temperature and pressure 
(which then is the appropriate form of the thermal equation of state), shows 
no dependence on pressure. This can be adapted to multicomponent mixtures, 
but such extensions will neither lead to constant total mass density nor to 
divergence free barycentric velocity fields; cf. HU. Concerning the mass density 
variations, a Boussinesq-type approximation seems possible but a rigorous and 
thermodynamically consistent modeling of such composition-induced volume 
effects is currently lacking. Such volume effects in the liquid case will be more 
severe in case of reactive mass transfer since the chemical reactions transform 
certain educts with their specific volume fraction into products which usually 
have different specific partial volumes. 

If chemical reactions take place in the liquid, for instance in bubble column 
processes, large sets of constituents are commonly involved. Despite the vol¬ 
ume effects mentioned above, this can lead to another significant complication 
occurring especially in aqueous solutions but also in other polarizable liquids: 
While the transfer components are uncharged, ionic species usually get involved 
during the course of chemical transformation. This couples the diffusive mass 
fluxed due to the appearance of Coulomb forces exerted by the electric field 
which is generated by the ions themselves. To account for these effects requires 
a coupling of the Maxwell-Stefan approach, say, to at least a quasi-static ver¬ 
sion of the Maxwell equations from electromagnetics, i.e. the Poisson equation 
to model the quasi-stationary electrical held generated by the ions. Since a 
detailed discussion of these model extensions is out of the scope of the present 
paper, we only refer to the recent publications HSU [3, but note that electrical 
effects can also play a role directly at huid interfaces. For the latter, see also 
the recent Ph.D. thesis m- 

The Outer Boundaries 

A sound model requires sensible conditions at the outer boundaries of the bal¬ 
anced region, which also includes a sound choice of appropriate outer boundaries 
itself. If possible, it is favorable to choose the outer boundaries to coincide with 
physical phase boundaries, like impermeable walls or a free liquid surface. In 
this case knowledge of the physics can be exploited to formulate reasonable 
boundary conditions in pretty much the same way as it has been done above 
for the fluid interface, i.e. by a combination of interfacial mass, momentum and 
energy balances with closure relations being motivated by the local entropy pro¬ 
duction rate. A severe complication for the mathematical description is the fact 
that the position of such physical outer boundaries is usually not fixed but will 
(slightly, say) vary due to the acting forces. In other words, the boundaries are 
usually pressure controlled, while almost all detailed mathematical descriptions 
employ volume control by hxing the spatial domain in which the considered 
processes take place. 
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In detailed or even direct numerical simulations one often cannot extend the 
domain up to physical outer boundaries, but has to introduce artificial bound¬ 
aries somewhere inside the bulk phases. At such non-physical boundaries the 
formulation of appropriate boundary conditions, so-called Artificial Boundary 
Conditions (ABCs), is a field in its own and an important topic in Computa¬ 
tional Fluid Dynamics. Let us note in passing that the use of perfect slip, i.e. 

v-n = 0 and (Sn)|| = 0, (73) 

models the boundary as impermeable which is in general wrong at an artificial 
inner surface inside a bulk phase. While quite some research has been devoted 
to the modeling of ABCs at outflow boundaries for single-phase flows (see, e.g., 
[23iini, the references given there and in [32]), much less is known for ABCs 
in two-phase flow situations. For rising or settling single fluid particles, an 
artificial boundary condition for lateral domain boundaries was derived in m 
which allows flow trajectories to exit and re-enter the computational domain. 
This allows for a reduction of the computational domain by about a factor of two 
in each lateral direction, thereby enabling a finer resolution near the interface, 
where it is crucial because of the HSNP mentioned above. 


Conclusion 

A multitude of physical phenomena can occur at fluid interfaces under strongly 
dynamical conditions. The potential of direct numerical simulations for gaining 
deep insights into such processes at fluid interfaces can only be fully exploited 
if the employed models provide a complete and sound description of the multi¬ 
scale multi-physics. For this purpose, the existing continuum mechanical models 
have to be critically assessed, improved and extended. The entropy principle 
provides a rational framework for such model extensions and its rigorous ex¬ 
ploitation leads to deep insights into, e.g., mass transfer across contaminated 
fluid interfaces. 
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Appendix: Interfacial Entropy Production 

We explain the main steps of a rigorous derivation of the entropy production 
in multi-component two-phase systems with non-vanishing interface mass den¬ 
sities. The starting point are the entropy balances in the bulk phases and on 
the interface, which read as 


dt{ps) + V • (psv -!-$)=( 


(74) 


in the bulk phases with s the specific entropy, $ the entropy flux and C the bulk 
entropy production, and 

dfip^s^) + Ve • -f -f |ps(v - v^) -f $1 • ns = (75) 

with the specific interface entropy s^, the interfacial entropy flux and the 
interfacial entropy production The second law of thermodynamics now 
states that 


> 0 and ^ 0 for any thermodynamic process. 


(76) 


The entropy densities are objective scalars which have to be modeled as material 
dependent quantities, where we consider the simplest set of primitive variables: 

ps = ps{pe,pi,...,pN) and p^s^ = p^s^{p^e^, pf,..., pf^). (77) 

The partial entropy derivatives have a clear physical meaning, which leads to 
the following definitions: 

dps Pi _ dps 1 dp^s^ pf __ dp^s^ 

“ T “ ^ ^ ~T^ 


1 

T 


'Y 


and 


1 

ys ■“ 


-• (78) 


dpi ' dp^e^ ’ ' dpf 

Using the chain rule, the balance equations (ITU) and (I75|) are used to express the 
entropy productions by means of the respective right-hand sides, where all time 
derivatives are eliminated by means of the corresponding balance equations. For 
this purpose the internal energy balance is required which reads 


N 


dt{pe) + V • (pev -I- q) = Vv : S -h ji • 


(79) 


i=l 


in the bulk phases with e the specific energy, q the heat flux and the individual 
body forces acting on the species Aj, and 


df{p e ) -I- Ve • (p e -|- q^) -|- {fnh + 


(v-v^)2 


- He 
N 


ns] 


(80) 


- [(v - v^)|| • (S™‘=nE)||] + [q ■ heI = Vev : -|- ■ bf 


i=l 


on the interface with the specific bulk enthalpy h = e -I- ^, the interface internal 
energy density p^e^, the interfacial heat flux q^ and the individual body forces 
bf. After lengthy but straightforward calculations, this yields 


N. N 

=V ■ ($ - ^ “ 7 f;ipe + P - psT piPi)V ■ V 


i=l 
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11 ^ y, 1 

+ q-V- + -Vv : S"- ■ (V^ - ^) - 
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(81) 
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for the bulk entropy production, where p is the thermodynamic pressure, Ra are 
the molar reaction rates of the a**' chemical bulk reaction and Aa the associated 
affinities defined as 


N 


Aa = 


(82) 


i=l 


with the stoichiometric coefficients z/f; see [0] for further details. For the inter¬ 
facial entropy production we obtain 


C^ = Ve-($^-^ + E^) 


(83) 
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The entropy fluxes are constitutive quantities which have to be modeled as 
objective vectors in such a way that the resulting entropy productions become 
sums of binary products, one for each dissipative mechanism. This leads to the 
constitutive relations 


N , -y N Y ,Y 

and = 

T T TS 




(84) 


for the entropy fluxes. Next, in order to satisfy the second law of thermody¬ 
namics, the term in front of the (surface) divergence of the (interface) velocity 
needs to vanish. This leads to the Gibbs-Duhem relations in the bulk and on 
the interface, i.e. 


N 


N 


pe+p- psT = y pipi and p^e^ - a - p^s^T^ = pf pf, (85) 


i=l 


i=l 


also called bulk, respectively interface Euler relations. Hence the reduced bulk 
entropy production 


11 ^ P, 1 

C = q • V- + -Vv : - y j, ■ (V^ - ^) - - y 

2=1 a—1 


( 86 ) 
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and, after specializing to vanishing interface viscosities, the reduced interface 
entropy production (1^^ result; see [3] for many more details on the entropy 
principle as the core part of continuum thermodynamics of chemically reacting 
fluid mixtures. For a recent extension to electrochemical interfaces see m- 
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